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A relevant problem regarding entanglement measures is the following: Given an arbitrary mixed 
state, how does a measure for multipartite entanglement change if general local operations are 
applied to the state? This question is nontrivial as the normalization of the states has to be taken 
into account. Here we answer it for pure-state entanglement measures which are invariant under 
determinant 1 local operations and homogeneous in the state coefficients, and their convex-roof 
extension which quantifies mixed-state entanglement. Our analysis allows to enlarge the set of 
mixed states for which these important measures can be calculated exactly. In particular, our 
results hint at a distinguished role of entanglement measures which have homogeneous degree 2 in 
the state coefficients. 



INTRODUCTION 

In recent years, there has been astounding experimen- 
tal progress in preparing and operating multi-qubit en- 
tanglement, e.g., with trapped ions [H0|, photons 0, Q|, 
and superconducting circuits 0, |(|. On the other hand, 
entanglement theory is lagging behind this advancement. 
To date, it is not possible to adequately characterize the 
entanglement of an experimentally prepared multipartite 
quantum state, that is, to quantify how much of a certain 
type of entanglement is contained in that state. 

One difficulty here is that it is not enough to just decide 
'how much entanglement is present'. Entanglement is re- 
garded as a resource for quantum information tasks, and 
different tasks require different resources, that is, inequiv- 
alent types of entanglement. The criterion whether or not 
two states can be considered equivalent with respect to 
their entanglement is whether they can be transformed 
into one another by means of stochastic local operations 
and classical communication (SLOCC) 0, S|- The con- 
sequence is that different equivalence classes of entangle- 
ment need to be distinguished pj-lllj. 

Entanglement witnesses provide a tool to distinguish 
SLOCC classes of multipartite states and to experimen- 
tally detect them [l2[. There are attempts to make wit- 
nesses more quantitative (l3l 14 1 which have not found 
widespread application in practice yet. 

So far there are only few practically relevant tools for 
measuring the entanglement of general (mixed) quantum 
states. Negativity-related entanglement monotones (TBI . 
[l6| represent one of them. They are applied to bipar- 
tite systems and it is not clear whether they can be used 
to distinguish different types of entanglement. There- 
fore we do not further investigate them here. The second 
well-established tool is Wootters' approach to compute 
mixed state entanglement for two qubits [l?} ■ It is based 
on the concurrence, a pure-state entanglement measure 
which is generalizable to pure states of larger systems 



than two qubits. The kind of pure-state entanglement 
measures the concurrence belongs to has several interest- 
ing properties (see below). In particular, entanglement 
measures of this kind can distinguish and quantify spe- 
cific types of multipartite entanglement. Recently there 
has been progress regarding the question how to obtain 
SLOCC classifications using them (i3 - l2fi |. An applica- 
tion of these quantities to mixed states in analogy to 
Wootters' method is therefore highly desirable. 

The concurrence is an example of a polynomial invari- 
ant, i.e., a polynomial function of the coefficients (and 
not their complex conjugates) of a quantum state, where 
the function is invariant under determinant 1 local oper- 
ations, that is under transformations 
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Here the operators Fk G GL(dk,C)® N represent general 
(invertible) local operations on the A-qudit Hilbert space 
H = <S>^ =1 C dh . Obviously, Ak have determinants equal 
to 1, i.e., A k e SL(d fc ,C). 

The peculiarity of polynomial invariants is that a 
power of them which has homogeneous degree 2 in the 
state coefficients is automatically an entanglement mono- 
tone [HI]. The generalizations of the concurrence, that is, 
polynomial invariants for pure states of three and more 
qubits are well-known |l9H-23j . 

Hence, pure-state entanglement measures are known 
in principle at least for multi-qubit systems. It is 
straightforward to generalize them to mixed states via 
the convex-roof extension (2tJ (for details, see below). 
The problem of quantitative entanglement theory is that 
it is not known how to compute the convex roof of an in- 
variant in general (or even a nontrivial lower bound of it) 
for any multipartite system, except for the concurrence 
of two qubits and its straightforward generalization to an 
even number of qubits, the N- concurrence (28l. [29|. 
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At least there exist a few exact solutions for three- 
qubit states of low rank and high symmetry (30l - l33t . To 
study how to extend these solutions to a larger family of 
states is part of the quest for a better understanding of 
the convex roof and multipartite entanglement in mixed 
states. An important question in this context is how an 
entanglement measure behaves under general local op- 
erations. While local unitary operations represent mere 
basis changes and do not change entanglement, general 
local operations require subsequent normalization of the 
state and thus a conversion of the amount of entangle- 
ment. 

In this work we investigate how multipartite entangle- 
ment needs to be rescaled if it is quantified by functions 
which are homogeneous and invariant under determinant 
1 SLOCC operations. It will turn out that in the case of 
mixed states this is easily done only for the convex roof 
of pure-state measures with homogeneous degree 2. That 
is, from the point of view of mixed-state entanglement 
rescaling one can conclude that there is a preferred ho- 
mogeneity degree of pure-state entanglement measures. 

The outline of this article is as follows. First we ana- 
lyze the rescaling properties of the convex-roof extension. 
Then we illustrate their application (and also their fail- 
ure for the 'wrong' homogeneity degree) by reconsidering 
some of the known exact solutions. Finally, we mention 
also some ways to extend the number of exactly solvable 
three-qubit problems. 



THE RESCALING METHOD 

This section provides the theoretical concepts on which 
this article is based. They are formulated for general 
homogeneous invariants, but we interpret them with re- 
gard to the physically relevant polynomial SL(dfc,C) 8Ar - 
invariants. To begin with, we introduce our notation and 
explain how a homogeneous invariant of pure states needs 
to be rescaled under a symmetry transformation if one ex- 
clusively focuses on normalized states. Then we explain 
the convex-roof extension. Finally, we state, derive, and 
discuss our main results. They describe the scaling of the 
convex-roof extension of a homogeneous invariant under 
a symmetry transformation for normalized mixed states. 



Rescaling of homogeneous invariants for normalized 
pure states 

We consider a real function r on some finite- 
dimensional Hilbert space H that is invariant with re- 
spect to some invertible linear operator L acting on 
%. In formulas, r(Lip) = t(i/j) for all ip G T~L. This 
implies t(-L _1 ^>) = T(tp). We assume further that r 
is homogeneous of degree r\. That is, for positive a, 
T(aijj) = a^r^ip). If we focus on normalized states <fr and 



4>' which are related by L according to <f>' = L(f)/\\L(j)\\, it 
follows that 



i<t>) 



(2) 



Note that t could be one of the polynomial invariants 
discussed in the introduction (or a function thereof) if 
H = 0f =1 C dk and L = A e SL(d k ,C)® N . Note further 
that Eq. (2| h as the same form for all rj such that there 
is no preferred homogeneity degree for the polynomial 
invariant. 



The convex roof extension 

Let H be a finite-dimensional Hilbert space and fl-u 
the convex set of (normalized) density operators acting 
on %. All p € fly, can be written as convex sums p = 
^iPi^i- That is, pi > 0, X^Pj = 1j an< ^ I s an extreme 
point of flfi . An extreme point tt of is given by ir = 
IV 1 ) (V' I / I Ml 2 ) where ip £ H (^ 0), and is usually referred 
to as pure state. A mixed state p has rank > 1 and its 
decomposition {pi,Tti} into pure states is not unique [34 1. 

A real continuous function / on H can be extended to 
flfi via convex-roof extension [27| : 



f(p) = ^m^Pifi^i)- 



(3) 



The minimum runs over all possible decompositions of p 
into pure states, and /(7Tj) = f{^i/\\ipi\\). A decompo- 
sition {pi,TTi} of p for which f(p) = YlPifi^i) * s called 
optimal. The minimal length of an optimal decomposi- 
tion (the minimal number of pure states in the decompo- 
sition) is < (rankp) 2 [27]]- Since /(/?) is an entanglement 
monotone on Q,n if / is an entanglement monotone on 
H [H, , the convex- roof extension is the standard way 
of applying a pure-state entanglement measure to mixed 
states. For furth er p roperties of the convex-roof exten- 



sion, see [271. l28l. |36|. However, as has been mentioned, 



finding the minimum in Eq. ([3|) is difficult in general. 



The convex-roof extension of homogeneous 
invariants 

Let us now investigate how the convex-roof extension 
of r transforms under the action of A, that is, how r(p) 
and r(p') are related if p' = LpL^ /tiLpL^ (p,p' € ^n)- 
Interestingly, it turns out that the convex-roof extension 
does not lead to a simple generalization of Eq. ([2]) for all 
rj. Rather, the convex-roof extension singles out invari- 
ants of homogeneous degree r\ = 2. Only in that case the 
convex-roof extension of r generally can be rescaled as 
for pure states, 



r(p') 



r(p) 
trLpLt 



(4) 
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This formula represents the case rj = 2 of a naive gener- 
alization of Eq. © to T(p') = T{p)/{trLpL^yi/ 2 which, 
however, is not correct in general. The main conclusions 
of this article are based on formula (@| . We will see that 
it can be understood as a consequence of the fact that 
L maps an optimal decomposition of p onto an optimal 
decomposition of p' if i] = 2. For r\ 7^ 2, the invariance of 
r under L still guarantees that r(p) and r(p') are either 
both zero or both nonzero if r(ip) > 0. These results 
give information about the convex roof of r for p' if it 
is known for p. Applying them to the polynomial invari- 
ants allows to extend (previous) results for the entangle- 
ment of a mixed state as measured by the polynomial 
invariants to all states that can be obtained from this 
state via general invertible local operations F (ILOs): If 
F = F 1 (g>...®F N e GL(d fc ,C) 8Ar , 



FpF^ ApAl 



trFpFt trApA^ ' 



(5) 



where A e SL(d kl C)® N is defined as in Eq. flTJ). An ex- 
ample that illustrates how Eq. Q can be utilized to cal- 
culate the three-tangle of mixed states will be presented 
in the next section. 

In order to see how Eq. (jU) is obtained suppose {pi, irA- 
and {qjiWj} are optimal decompositions of p and p' of 
lengths m and n, respectively. In other words, r(p) = 

YnLiPi T ( 7T i) and r G°') = YJ]=il3 T ( w j)- We express p' 
in terms of the {pi,TTi} and p in terms of the {qj,zaj}, 



i=i j=i 



P 



(6) 
(7) 



where 



p\ = piT-hrLntf, q'j = q J Ttx[L- 1 w j (L- 1 )t], (8) 
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We have used the abbreviation T = tvLpLK It is easy to 
see that p[ = q'j = 1 and that 71^ and hj^ are pure 

states. Since LwiL^(= LiTiirjL^) is a positive operator 
and L is invertible, trL7riL^ > 0. Consequently, T > 0. 
A similar argument ensures that tr[L" 1 n7 i ( J L- 1 )t] > 0. 
Starting from Eqs. © and ([7]), one can estimate 



Tip) < Tj^^L-^iL-y^T^ 



.))■ 



r(p') < T-^pilixbKiL^Tfa). 



(10) 



(11) 



If r/ = 2, Eq. (0| follows. Moreover, one sees that {p^ir^} 
is an optimal decomposition of p' because ^iPiTfjr'A = 



^2iPiT{iTi)/T = r(p)/T. Further, {q'^w'j} is an optimal 
decomposition of p. Hence, an optimal decomposition 
{pi, TTi} {{qj,vjj}) of p (//) is mapped by L (L -1 ) onto an 
optimal decomposition of p' (p) which is given by {p' t , ir^} 
{{q'j^'j}) in Eqs. © and ©. 

If r\ 7^ 2, but t > 0, one can still conclude from Eqs. 
(fTOf and (fTTj) that either both r(p) and r(p') are zero 
or both are nonzero. The assumption r(p) = implies 
t(tt 1 ) = for all m and, due to Eq. ([IT]), r(p') = 0. On 
the other hand, the assumption r(p) > implies that 
r(vjj) > for at least one Wj because of Eq. (|TD|) . The 
converse follows analogously. 

Let us briefly discuss these results with regard to the 
polynomial invariants and their application as entangle- 
ment measures. As already indicated, our results, in par- 
ticular Eq. ((3]), provide a useful tool for learning about 
the convex-roof extension of a polynomial invariant for a 
mixed state p' if it is already known for a state p that 
can be transformed into p' via SLOCC. However, we have 
also gained another important insight about the convex- 
roof extension of polynomial invariants: It has peculiar 
properties for polynomial entanglement measures of ho- 
mogeneity degree 2 (or, of linear degree if the pure-state 
measures are regarded as functions of rank-one density 
operators) and therefore indicates a special significance 
of the latter. Only for entanglement measures of this 
type, the invariance property for normalized pure-states 
generalizes to normalized mixed states and the entangle- 
ment both of pure and mixed SLOCC-equivalent states 
is related by a simple analytical formula, Eq. (j4|). 



EXAMPLES 

Now we consider several examples to illustrate how 
the rescaling method from Section n can be applied to 
generalize existing results for mixed-state entanglement. 
We also will demonstrate what happens if one tries to 
apply it to monotones of degree 7^ 2. 

The examples we consider are mixed states of three 
qubits. For three-qubit systems, there is only one poly- 
nomial SL(2, C)® 3 invariant, the three-tangle [ljj. For 
W) = Elj-,fe=o^t \ijk) it is given by 

r 3 (V>) = 4|di-2(fe + 4d 3 |, 

d i = V'oooV'iii + V'ooiV'iio + V'oioV'ioi + V'on^ioo 
d-2 = V'oooV'ooiV'iioV'm + V'oooV'oioV'ioiV'iii + 
+^000^011^100^111 + V'ooiV'oioV'ioiV'iio + 

+V'00lV ; 01lV'l00''/'ll0 + ^OlO'0OllV'lOO''/'lOl 

d3 = ^oooV'iioV'101^011 + ^iooV'oioV'001^111 (12) 

and distinguishes Greenberger-Horne-Zeilinger (GHZ) 
entanglement from W-type entanglement. For mixed 
states, it is defined by convex roof extension. 

The three-tangle has homogeneous degree 4. However 
to apply the rescaling method Eq. (01 as described above, 
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degree 2 is required. Therefore on pure states we have to 
use the square root of the three-tangle, r(tp) = \/t^(iI>) 
which is an entanglement monotone as well |37|. For 
mixed states, we again define r(p) through the convex 
roof extension. Note that this is not the same as taking 
the square root of Ts(p). However r(p) = iff r 3 (p) = 0, 
therefore both are equally suitable to distinguish GHZ- 
type and VF-type entanglement. 



Mixtures of generalized GHZ and W states 



We consider r(p) for rank-2 mixtures 

p(p) = V |gGHZ) (gGHZ| + (1 - v) W) {gW\ (13) 

of a generalized GHZ and generalized W state (both nor- 
malized) 



IgGHZ) = a|000) + 6|lll) 
\gW) = c|100> + d|010> 



|001) 



(14) 



The three-tangle of those states has been calculated in 
[3ll | . We will be able to reuse pure-state results from 
that paper. 

We begin by solving the problem using the character- 
istic curve method [38( and then demonstrate how those 
solutions can be mapped to each other by means of the 
rescaling method. 

The characteristic curve is defined as the minimum 
tangle of the states at the same "height" in the Bloch 
sphere as the corresponding mixed state, 



f(p) = min(r(v^ |gGHZ) + y/T^j* \gW))). (15) 

4> 

Since we are dealing with pure states at this point, and 
the square root is monotonic, we can reuse the results of 
Ref. [31| and get 



T(p)=2\ab\J p 2-\ s \^(T^pf 

where we have used the definition 

icdf 



(16) 



(17) 



The convex characteristic curve T*(p) is the function 
convex hull of f(p), i. e. the largest convex function which 
is nowhere larger than f(p). As was shown in Ref. [HI, 
the convex characteristic curve is always a lower limit to 
r(p). Using the results in [3~0 | we find that f(p) has zeros 
at p = and at p = po with 



S3 



Po = 



1 



(18) 



It is easy to check that for p > po, f(p) is concave (note 
that this is different from the three-tangle itself where 



fs(p) is convex right above po for s < 2\/2). Therefore in 
this region, the convex characteristic curve is a straight 
line. We thus get 



T*(p) 







for < p < po 



2\ab\[p-(l-p)\s\3) forp <p<l 



(19) 



In the current case, the convex characteristic curve ac- 
tually gives the correct value of r{p). This can be seen 
from the symmetry of the problem: Applying the local 
unitary transformation 



U 3 



cxp 



exp 



r (2) 



cxp 



r(3) 



(20) 



changes neither the GHZ nor the W state (and thus also 
not their mixtures), nor, of course, the three-tangle. But 
it changes the relative phase between the GHZ and the W 
state in the superpositions of Eq. ([T5]) , thus rotating the 
Bloch sphere spanned by these superpositions by 2tt/3. 
Thus each state gives rise under this symmetry U3 to 
three locally SU equivalent states whose equal mixture 
gives exactly p(p) (Eq. (fT3"|) ). This is especially true for 
the minimum r state defining the characteristic curve. 
Thus, due to the symmetry the characteristic curve gives 
an upper limit to r. 

Since the convex characteristic curve always provides a 
lower limit, wherever both agree the unconvexified curve 
gives the correct tangle, and therefore in those points the 
above decomposition is an optimal one. 

While the above argument applies only to the points 
p = 0, p = po and p = 1, it is obvious that for the lin- 
ear sections of the convex characteristic curves, the value 
r*(p) can be achieved by appropriately combining the 
optimal decompositions of the corresponding end points 
which proves that indeed T(p(p)) = r*(p). 

After deriving the solution of the problem, we now 
show that different states (TIB")) are SLOCC equivalent 
and, hence, that the solutions r for different values of 
a, . . . , / can be obtained from each other via the rescal- 
ing method. To this end, we apply a diagonal invertible 
local operation of determinant 1, 



A 



a 
l/a 



(3 
1/(3 



7 

1/7 



(21) 



which transforms |gGHZ) to M |GHZ) and \gW) to 
A/2 I W) ■ An easy calculation shows that this is achieved 
for 



adf 







acf 



1 



bp 

acd 



(22) 



This transforms |gGHZ) into v^lGHZ) and \gW) to 

1 

v/3 [ \l\cdf) 3 \W) = i/| • V2ab ■ (2s) 3 \W) and there- 
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fore p into 

A P A* =2\ab\ fp|GHZ)(GHZ| + 3 g ^ |2«|* \W)(W\j 

(23) 

Since the trace of this operator is 

T = tr(ApAl) = 2 |o6| (p + \2s\ 2/3 \ , (24) 



we have to divide p' by this expression and find 



T 



p'\G'BZ){GSZ\ + [l-p')\W){W\ (25) 



with 



or solving for p: 



P = 



,+ Ml_£)|2 S |2/3 



3p'\2s\ 2 / 3 
3p'|2s| 2 / 3 + 8(1 - p') 



(26) 



(27) 



Note that this is a monotonously increasing function of 
P'- 

Now we can apply the rescaling method described in 
the previous section to this problem. By using Eq. (|4]) 
we obtain 



t(p(j>)) = T T(p'(p')) 



(28) 



where T is taken from Eq. (|2"4"]) and p' from (|2l)]) . 

For the standard GHZ/IT problem, that is a = b and 
c = d = /, we have s = 2 7 / 2 /3 3/2 and 2|a6| = 1, and 
therefore 



r( P '(p')) 



I for < p' < p' Q = 

\p'-(l-p>)^ {ovp' <p'<l 



2 7/3 

3+2 7 /3 



(29) 

By inserting all values one indeed recovers Eq. (fT9|) with 
the correct value of po, Eq. (fT8")) , see Fig. [T] This means 
in turn that, by virtue of the rescaling method, the knowl- 
edge of the solution of the standard problem with sym- 
metric GHZ and W states in Eq. ([29| is sufficient to ob- 
tain the solution of the more general states in Eq. (fT3| . 
This is achieved simply by applying ILOs according to 
Eq. (f2"Tj) and appropriate rescaling following Eq. JU) . 

Note that while we have considered only diagonal 
transformations the method is of course valid for general 
invertible local operations. Those break the symmetry 
noted above, and therefore lead to cases which cannot be 
solved using the convex characteristic curve method. For 
those problems, no general method is yet known. 



Using the rescaling method for degrees other than 2 

As mentioned before, the rescaling method is not fully 
applicable to entanglement measures with homogeneous 



GHZ 



gGHZ 




A 




FIG. 1: Illustration of the rescaling method for the square 
root of the three-tangle r of mixtures of standard (left) and 
generalized (right) GHZ and W states. The two orthogonal 
eigenstates (g)GHZ and (g)W are located at the north and 
south poles of the Bloch spheres, respectively, and the mixing 
parameters p' , p vary along the lines connecting the two poles. 
The symmetric GHZ/W problem (left) can be transformed 
into the generalized GHZ/IT problem (right) and vice versa 
by an invertible local operation. The corners of the tetrahe- 
dra show the pure states of the optimal decomposition of the 
states inside those tetrahedra. In the lower (blue) tetrahe- 
dron, t — 0, while in the upper (green) tetrahedron, it grows 
linearly with p' (p). If the convex-roof extension is known 
for one of the problems, it can be calculated for the other by 
virtue of the rescaling method. 



degree rj ^ 2. However, by considering the three-tangle 
(j) = 4) of GHZ/W mixtures we will see that if the prob- 
lem has certain symmetries, the rescaling method can 
still be useful for calculating entanglement measures with 

The three-tangle of the mixtures in Eq. (fT3")) has al- 
ready been calculated in [3l[, and for the special case of 
mixing standard GHZ and standard W in [3(| (see also 
Fig. [2|) . For the mixtures of standard GHZ and standard 
W, the three-tangle has the form 



rsip') = 



-^P'(l-P'f 
(1-p' 




0<p'<Po 
<p'< Pi 
< p' < 1 



Pi 

(30) 

where p' = 3 ^_ 27 / 3 and p\ = \ + We now derive 

the results of a simple-minded application of the rescaling 
method to the 3-tangle which is based on the naive gen- 
eralization of Eq. © (below Eq. (TJ}) with r/ = 4. That 
is, we calculate p' from p using Eq. (|26jl , apply the func- 
tion above, and then multiply with T 2 from Eq. (|24|) . It 
will turn out that this procedure indeed leads to incorrect 
results in general. This is also suggested by Fig. [21 Opti- 
mal decompositions of states which can be transformed 
into one another by means of operations Eq. f|2 1 1) can 
have different lengths. Hence, the ILOs in general do 
not map optimal decompositions onto optimal decom- 
positions, and one would therefore not expect that the 
three-tangle of a mixed state can be traced back to the 
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(known) three-tangle of an SLOCC-equivalent state via 
the rescaling method. However, we will also see that in 
some cases the correct three-tangle of mixtures of gener- 
alized GHZ and W states can be obtained from Eqs. (|5D|) 
by means of this method due to the particular symme- 
tries of this problem. 

To identify the range of p where the three-tangle van- 
ishes we insert p'n into Eq. (|27| and get Eq. (fl8| in agree- 
ment with Ref. [31J. Indeed this was to be expected, be- 
cause the method allows to reliably distinguish between 
zero and non-zero values of the polynomial even for ho- 
mogeneity degrees other than 2. 

Next, we calculate the tangle for p < p < p\ using 
the method above, where Pi is the p value beyond that 
r 3 (p(p)) increases linearly [31|. From 



T 3 ( P (p))=T 2 T 3 (p'(p>)) (31) 



r z {p)^A\ab\ 2 (p 2 ~\sWp(l-pf) (32) 



we get 



which again agrees exactly with the result in Ref. [31|. 
This comes somewhat as a surprise, because the method 
should not be applicable here. However, it is easily ex- 
plained by considering the symmetry of the problem. 

Let us suppose the special case that both {pi^i} and 
{pi, AiTiAT /T} are optimal decompositions of p and p' in 
normalized pure states. This implies ^(AiTiA^) = T for 
all i. One easily sees from Eqs. (| 10[) and (|TT]) that under 
these conditions r(p) = T''/ 2 t(//) if r is an arbitrary 
homogeneous invariant, or, in case of the three-tangle, 



gGHZ 



GHZ 



gGHZ 




FIG. 2: Illustration of the three-tangle tz of mixtures of stan- 
dard GHZ and W states and generalized GHZ and W states. 
As in Fig. [T] the (g)GHZ and (g)W states are at the north 
and south poles of the Bloch spheres, respectively, and the 
corners of the tetrahedra and triangles give the elements of 
the optimal decompositions. In each red triangle, tz = const. 
In the blue tetrahedra, tz — 0, and in the green tetrahe- 
dra tz depends linearly on the mixing parameter. Here, the 
sphere in the center represents the standard GHZ/IF case. 
The transformation to the generalized problem leads to two 
different cases: For s < 2\/2 (left sphere), the structure of 
the symmetric case is preserved. However, for s > 2\/2, the 
two tetrahedra sit directly on top of each other. As states 
equivalent under ILOs can have optimal decompositions for 
tz with different lengths, the figure indicates that the rescal- 
ing method cannot be applied to calculate tz- 



Tj,(p) = T 2 Tj,(p'). The specific mixed states and ILOs of 
our example satisfy these conditions for po < p < p\ , as 
we now explain. 

In this region characteristic curve and convex charac- 
teristic curve agree with each other. As described above, 
this is because there is an optimal decomposition which 
consists purely of states related by the local symmetry 
(|2T)|) which rotates the Bloch sphere by 27r/3. 

Arbitrary rotations of the Bloch sphere about the 
GHZ/IT axis cannot be achieved by local operations, 
however, they can be realized by non-local diagonal uni- 
tary transformations, e.g., U(tp) = diag(l, e lv , e lv , e~ lv , 
e iv , e" iv , e~ iv , 1). Note that the local operations (J2UJ) are 
a subset of these. Now the ILOs (f2"Tj) are also diago- 
nal, and therefore commute with those transformations. 
This means especially that if two vectors are equivalent 
under this symmetry U(tp), the transformed vectors are 
also equivalent under this symmetry, which means that 
horizontal planes in the Bloch sphere are just moved ver- 
tically by the transformation. Moreover, it is easy to 
verify that if two pure states can be transformed into 
each other by U(ip), i.e., share the same latitude, also 
their norms are multiplied by the same factor. This 
implies that for those states the transformed state has 
minimal three-tangle iff the original state had minimal 
three-tangle. Of course, all convex combinations of those 
states are also multiplied with the same factor. Espe- 
cially, the characteristic curve of the original problem is 
mapped onto the characteristic curve of the transformed 
problem. 

However, to get the correct three-tangle, one needs the 
convex characteristic curve, which does not agree with 
the characteristic curve above p\. Now it is not to ex- 
pect that the transform of the latter gives the same result 
as the convexification of the transformed curve, and that 
optimal decompositions are mapped to optimal decom- 
positions in general. And indeed, if one tries a straight- 
forward calculation of p\ using Eq. (|27|1 one obtains 



Eq. j2Tl 

Pi 



24 |2s|f 
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24 |2s|f 



(33) 



which not only does not agree with the correct p\ from 
(3l| . but even has the wrong behavior. While the correct 
Pi decreases from 1 to 1/2 on increasing s (actually, only 
until it hits pa), p^ q ' instead increases from to 1. 

As the unconvexified characteristic curve is mapped 
onto the unconvexified characteristic curve of the trans- 
formed problem, the correct three-tangle can be obtained 
by using the former to calculate the convex characteristic 
curve directly. However, we emphasize again that this is 
only due to the high symmetry of the problem and the 
fact that the symmetry operations commute with the di- 
agonal ILOs. 

As in the case of degree 2 (Section p), one might con- 
sider to further extend the range of solutions derived from 
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the standard GHZ/W case by applying arbitrary (i.e., 
not only diagonal ones) local SL transformations. How- 
ever, general transformations modify the unitary sym- 
metry of the problem into one which is non-unitary, 
AU(ip)A~ 1 . Such a transformation in general does not 
map the characteristic curve (and the optimal decom- 
positions) of the standard GKZ/W mixtures onto the 
characteristic curve (and the optimal decompositions) of 
the transformed mixed states. Therefore, in general the 
rescaling method is not useful for finding new solutions 
for the convex roofs of entanglement measures with ho- 
mogeneity degree r/ ^ 2. However, note that the zero 
polytope of the original problem will be mapped exactly 
onto the zero polytope of the transformed problem, irre- 
spective of rj (see Sec. J. 



Using the rescaling method for mixed states of rank 
greater than 2 

If the degree of homogeneity of the polynomial invari- 
ant is 2, the rescaling method is fully applicable. We 
show now by discussing further examples that it can 
be also successfully employed for calculating the entan- 
glement of mixed states with rank > 2. Consider, for 
instance, rank-3 mixtures of |GHZ), \W) and \W} = 
(| 011) + 1 101) + 1 110))/ V3. The three-tangle for such mix- 
tures has been found by Jung et al. in Ref. [32| ■ We illus- 
trate the application of the rescaling method for states 
of rank > 2 by briefly sketching how the square root of 
the three-tangle can be obtained for those states and all 
SLOCC e quiv alent states. Adopting the parametrization 
from Ref. [32j, we have 

p(p,q) = p\GUZ)(GUZ\+q\W)(W\ + (l-p-q) \W)(W\ 

(34) 

The analog of the characteristic curve is now the two- 
dimensional "characteristic surface" 



because the surface is not convex 
Eq. 



are the ones given in 



r(p, q) 



min r[v^|GHZ) 



-Vqe^ \W) - 



qe 1 ^ 2 \W) (35) 



where r = \Jtz. We have a gain adopted the parametriza- 
tion of the states used in [32j. As before, the problem 
has the symmetry (|20|) . and an equal mixture of a set 
of states equivalent to a given state under this symme- 
try states gives p(p,q). It is easy to see that r takes its 
minimal value for 4>\ = 4>2 = 0, which results in 



t(p, q) 



p2 _ 4^ 



(36) 

where we have defined r = 1 — p — q. The decomposi- 
tions of p(j>, q) leading to those values (note that in gen- 
eral they do not represent the optimal decompositions 



11) of [32 1 . 

When plotting this function (see Fig. [3]), one easily sees 
that there is a curved line of zeros. All states between 
that line and the connecting line between \W) and \W) 
are W-type states with r = T3 = 0. One also sees in 
the graph in Fig. [3] that beyond this line of zeros, one 
can convexify the surface by just connecting the points 
of the zero line with the GHZ state using straight lines. 
Calculating the exact values lies beyond the scope of the 
present paper as our aim merely is to illustrate the po- 
tential of the application of the rescaling method. We 
note, however, that several points on the line of zeros are 
given in table I of Ref. [32] . Once the solution for the 
convex roof of r = y/r^ is obtained, it is straightforward 
to apply the criterion |HJ|. This allows to calculate the 
square root of the three-tangle for a larger set of states. 

We conclude this section by noting that also in the 
case of these rank-3 mixtures even the three-tangle T3 
(degree 4) could be extended by a modified application 
of rescaling. This is (in complete analogy with the previ- 
ous section) because the states in (|35|) whose T3 is mini- 
mized for the characteristic surface are equivalent under 
a diagonal (non-local) unitary symmetry. Therefore, by 
considering only diagonal ILOs, the characteristic surface 
of the three-tangle (without square root) is again mapped 
to the characteristic surface of the transformed problem. 
Thus in analogy to the preceding section we can calculate 
the three-tangle "after ILOs" also for this case from the 
three-tangle of the original problem, by convexifying the 
transformed surface. 

Finally, we remark that there are cases of even higher 




0.5 T 



0.0 



FIG. 3: Characteristic surface f(p,q) (Eq. (36)) for mixtures 
of standard GHZ, standard W and flipped W state. As de- 
scribed in the text, due to the symmetry of the problem the 
convex hull of this function gives the correct value of r(p). 
The zeros of r are found for the W state (p = 0, q = 1), the 
flipped W state (p = 0, q = 0) and for a continuous curve 
starting at p — q — p' (see Eq. (27)) and ending at p — p' , 
q — 0, where the characteristic surface has a kink. The con- 
vex hull of this function is zero between the line connecting 
Wyflipped W and the curved line of zeros described before. 
From the graph it is also easy to see that between that curve 
and the GHZ state, the convex hull consists of the straight 
line segments which connect the individual points of the curve 
with the point corresponding to the GHZ state. 



s 



rank where the entanglement of classes of SLOCC-equi- 
valent states can be obtained via the rescaling method. 
The three-tangle of certain mixtures of three-qubit GHZ 
states of high rank has been calculated by Shu- Juan et al. 
[33l | . It should be straightforward to carry out also those 
calculations for the square root of the three-tangle. Then 
the rescaling method can be used to extend the scope of 
those results. Since Shu- Juan et al. used the convex 
characteristic curve method, it should also be possible 
to calculate, for a suitable restriction of the local trans- 
formations, the characteristic surface of the transformed 
problem for the three-tangle. Its convexification should 
yield the three-tangle of the transformed problem. 



CONCLUSION 

We have developed a rescaling method to calculate 
polynomial-based entanglement monotones of homoge- 
neous degree 2 (in the coefficients of a pure quantum 
state) in mixed states. The method is based on trans- 
forming the mixed state under consideration into one of 
an already solved problem by using invertible local opera- 
tions. We have further demonstrated that rescaling is not 
generally applicable to functions of homogeneous degrees 
other than 2. Therefore, from the point of view of the 
convex-roof construction, there is a clear preference for 
entanglement monotones of homogeneous degree 2 (and 
correspondingly, of homogeneous degree 1 in the coeffi- 
cients of the density matrix). This is very much in line 
with some conclusions by Verstraete et al (l8| . 

The application of the rescaling method has been illus- 
trated by means of several examples. We have found that 
the presence of a certain symmetry in some cases even 
allows to calculate the convex roof of entanglement mea- 
sures with homogeneous degrees other than two. More- 
over, the question whether the monotone is zero or non- 
zero, which is important for SLOCC classification, can 
be reliably answered with this method independent of 
the homogeneous degree. 
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